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Abstract 

The reorientation transition of the magnetization of ferromagnetic 
films is studied on a microscopic basis within Heisenberg spin mod- 
els. Analytic expressions for the temperature dependent anisotropy 
are derived from which it is seen that the reduced magnetization in 
the film surface at finite temperatures plays a crucial role for this tran- 
sition. Detailed phase diagrams in the temperature-thickness plane are 
calculated. 

The direction of the magnetization of thin ferromagnetic films depends on 
various anisotropic energy contributions like surface anisotropy fields, dipole 
interaction, and eventually anisotropy fields in the inner layers. These com- 
peting effects may lead to a film thickness and temperature driven spin re- 
orientation transition (SRT) from an out-of-plane ordered state at low tem- 
peratures to an in-plane ordered state at higher temperatures at appropriate 
chosen film thicknesses. Experimentally, this transition has been studied in 
detail for various ultra-thin magnetic films [||, ||, Recently, it was found 
by Farle et al. [|] that ultra-thin Ni- films grown on Cu(OOl) show an opposite 
behavior: the magnetization is oriented in-plane for low temperatures and 
perpendicular at high temperatures. 

It has been shown by various authors that the mechanism responsible for 
the temperature driven transition can be understood within the framework of 
statistical spin models. While Moschel et al. || showed that the temperature 
dependence of the SRT is well described qualitatively within a quantum 
mechanical mean field approach, most other authors focused on classical 
spin models. Results obtained by extended Monte Carlo simulations || [7]] 
as well as mean field calculations of classical Heisenberg models (|, || |10| , 
although different in detail, agree in the sense that a temperature driven 



Journal of Magn. and Magn. Mater., accepted for publication (1998) 



2 



SRT is obtained. Very recently, the perturbative approach |10[ has been 
extended to itinerant systems [|ll|]. Note, however, that due to an unsufficient 
treatment of the spin-orbit induced anisotropy in this work no SRT is possible 
for mono- and bilayers in disagreement with the above mentioned theories 
and with experiments. 

While the SRT from a out-of-plane state at low temperatures to an in- 
plane state at high temperature is found to be due to a competition of a 
positive surface anisotropy and the dipole interaction the interesting new 
result for ultra-thin Ni-films is argued [|J to have its origin in a stress-induced 
uniaxial anisotropy energy in the inner layers with its easy axis perpendicular 
to the film. This anisotropy is in competition with the dipole interaction and 



a negative surface anisotropy. Theory confirms this picture |9|, 10 . 

In this paper a systematic investigation of the parameter and thickness 
dependence of the reorientation transition and in particular a calculation of 
the corresponding phase diagrams is presented. The calculations are done in 
the framework of a classical ferromagnetic Heisenberg model consisting of L 
two-dimensional layers on a simple cubic (001) or face centered cubic (001) 
lattice. The Hamiltonian reads 

ft = - « s i ■ s i ~ £ D M ) 2 + 7f Yj r i? s i • s i ~ 3r ij 5 ( s i ■ r *i)( r *i ■ s i)> (!) 

Z (ij) i Z V 

where s, = (sf, sf , sf) are spin vectors of unit length at position i-j = (rf , rf , rf) 
in layer \ and = — r^. J is the nearest-neighbor exchange coupling 
constant, D\ is the uniaxial anisotropy which depends on the layer index 
A = 1...L, and lj = /x /x 2 / Ana 3 is the strength of the long range dipole 
interaction on a lattice with nearest-neighbor distance a (/xo is the magnetic 
permeability and \i is the effective magnetic moment of one spin). Note that 
only second order uniaxial anisotropics D\ enter the Hamiltonian ([!]). In our 
calculations we will restrict ourself to the case that all anisotropies are equal 
to D v except at one surface, where the value is D s , as this scenario is suffi- 
cient for explaining the basic physics of the temperature driven reorientation 
transition. 

The Hamiltonian ([J) is handled in a molecular-field approximation, the 



reader is referred to |10[ for further details. In the following we assume trans- 
lational invariance within the layers and furthermore that the magnetization 
is homogeneous inside the film and only deviates at the surfaces, which is 
justified by numerical simulations [p]0| . Therefore we can set (sj) = (s^) if Sj 
is a spin in a volume layer (A = v) or in a surface layer (A = s). Now we 
must distinguish two cases: 

• For L <= 2 we get a system with only one (surface) layer with the 
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effective interactions x ss = x$ + 2x\{l — Lj 1 ), where x is either the 
coordination number, z, or the dipole sum, respectively. In this 
case we have L s = L and L v = 0. The effective uniaxial anisotropy is 
D s = D s + (L s — 1)D V and the Curie temperature of the unperturbated 
system is T C (L)/J = z ss /3. 

• For L > 2 we get a system with two layers with the effective interactions 
x ss = x , x vs = xi, x vv = xq + 2xi(l - L~ x ), and x sv = 2xiL^ 1 . Here 
we have L s = 2 and L v = L — 2. Note that the coupling between the 
volume spin and the surface spin is asymmetric. The effective uniaxial 
anisotropics are D s = D s + D v and D v = L v D Vl respectively. The 
critical temperature of the unperturbated system is given by 

In this approach we are not restricted to films with integer values of the 
thickness. The number of next neighbors z§ and the constants containing 
the effective dipole interaction between the layers for different lattice types 
and orientations are listed in table [l|. 



lattice 


zo 




ZS>1 


$ 




$5>1 


sc(OOl) 


4 


1 





9.034 


-0.327 




fee (001) 


4 


4 





9.034 


1.429 


( e -V2*(«-l)) 



Table 1: Numerical constants z$ and for sc(001) and fcc(001) lattices. 

A perturbation theory to the mean field Hamiltonian can be applied con- 
sidering the terms proportional to u and D\ as small perturbations which is 
justified in view of the smallness of these parameters || [IIJ. In this frame- 
work we derived an analytical expression for the total anisotropy K(T) of the 
system which is defined as the difference of the free energies of the in-plane 
state and the out-of-plane state. K{T) depends on temperature only through 
the absolute value of the subsystem magnetizations m\{T) calculated with 
the unperturbed Hamiltonian and is given by 

^) = £^(l-3T!^) ~Em^Lvm y (r). (2) 

hA = JJ2\> Z\\i th\i is the mean field in subsystem A. If we neglect the narrow 
canted phase, the reorientation temperature T r is given by the condition 
K(T r ) = 0. 
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For the homogeneous case [L <— 2) the magnetization is approximately 
given by 

2,^ l-T/T e 

m 2 (T) « 

y ' 1 - 2T/5T C 

This approximation deviates less than 1% from the values obtained numeri- 
cally. From (d) we then get 



K(T) « ( 1 - £ 



4(1 - 2T/5T C 



from which the reorientation temperature is determined as 

T r _ 5 15w$ L + 2$ 1 (L- 1) 
T~ c ~2~ 8 A, + (L — 1) A, 

The phase diagram of the fcc(OOl) system calculated within this approach is 
shown in figure [I]. One can distinguish three different reorientation schemes: 

• For uniaxial volume anisotropics 

— < -($ + 2$i) -8.92 (fee) 
uj 4 

we find a normal SRT from perpendicular to in-plane magnetization 
with increasing temperature and film thickness (solid lines in figure |l]). 

• If D v is large, we find a reversed SRT from in-plane to perpendicular 
magnetization with increasing temperature and film thickness (dashed 
lines in figure |l|) . 

• In the intermediate range we obtain a third type of SRT where the mag- 
netization switches from perpendicular to in-plane direction with in- 
creasing temperature, but with decreasing film thickness, (dash-dotted 
lines in figure |1|) . 

This work was supported by the Deutsche Forschungsgemeinschaft through 
Sonderforschungsbereich 166. One of us (A.H.) would like to thank M. Farle 
for fruitful discussions. 
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Film thickness [ML] 

Figure 1: Phase diagrams of the fee system for several values of D v /uj 
(Numbers). D s /uj is chosen so that in the ground state the SRT occurs 
at L r (T = 0) = 6, i.e. D s = 51.37^ — 5-0^. The thick solid line represents the 
Curie temperature T c of the film, the solid lines are normal (Fe-type) SRTs, 
while the dashed lines are reversed (Ni-type) SRTs. The dash-dotted lines 
represent a third form of SRT, where the reorientation temperature T r grows 
with film thickness and does not necessarily reach T c for finite values of L. 
At the thin dotted line D v /u — > ±oo, i.e. here the SRT occurs for uj = 0. 



